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■ Using the Landshoff-Nachtmann model of the pomeron, rapidity and 

■ transverse momentum distributions of gluon jets produced by double 
pomeron exchange in pp (pp) collisions are calculated. The comparison 
with the CDF Run I results on the central inclusive dijets production 
cross sections is performed. We find the model to give correct order of 



magnitude for the measured cross sections. 
PACS numbers: 13.87.Ce, 13.85.-t, 12.38. Qk 



1 Introduction 



The subject of Higgs boson production by double pomeron exchange (DPE) has 
^ \ drawn noticeable interest in recent years [1—11] . 



The lack of solid QCD framework for diffraction makes the determination diffi- 
cult. Despite the distinct progress in the recent years the serious uncertainties are 
still present. 

One way to reduce these uncertainties is to study other double pomeron exchange 
processes and compare them with existing data. A particularly illuminating process 
is the DPE production of two jets (dijets). Such a process was originally discussed 
at the Born level in [12]. Later the dijets production was studied in [5, 13] and in 
[8-11,14]. 

Recently, using the Landshoff-Nachtmann model of the pomeron, the cross- 
section for gluon jets production was calculated [15]. The obtained results together 
with those for quark-antiquark jets calculated some time ego [16, 17] give the full 
cross-section for dijet production in double pomeron exchange reactions 1 . 

1 It is generally accepted that DPE dijets mainly consist of gluon jets. For that reason we will 
use gluon jets and dijets interchangeably. 
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In the present paper we calculate the differential cross-section for gluon jets 
production hoping that this will allow a more precise comparison of the model with 
experiment. 

The comparison with the CDF Run I results on the central inclusive dijets pro- 
duction cross sections [18] is also performed. 




Figure 1: Production of two gluons by double pomeron exchange. 



2 Matrix element and notation 

In the Landshoff-Nachtmann model [19] the pomeron is approximated by an ex- 
change of two non-perturbative gluons coupled to one of the quarks of the colliding 
hadrons. 

The matrix element for gluon jets by double pomeron exchange in such model is 
given by the sum of the three diagrams shown in Fig. 2 (plus analogous emissions 
from the second gluon) 




Figure 2: Three diagrams contributing to the amplitude of the process of gluon pair 
production by double pomeron exchange. The dashed lines represent the exchange 
of the non-perturbative gluons. 

where the inner quark lines are put on shell. In this model, it was shown that the 
square of the matrix element (averaged and summed over spins and polarizations) 
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for the diffractive production of two gluons is of the form [15]: 



\M pp \ 2 = %l\M qq \ 2 [F{t l )F{t 2 )} 2 , (1) 
where |M gg | 2 is the dijet production amplitude squared for colliding quarks 2,3 : 

= C ~, £} ,2 %~ Mtl) ^l Mt2) exp (2(3 ( tl + t 2 )) R 2 . (2) 

Transverse momenta of the produced gluons are denoted by u 1 and u 2 . The constants 
C and R are defined later, a (t) — 1 + e + a't is the pomeron Regge trajectory with 
e « 0.08, a' = 0.25 GeV~ 2 (t 1: t 2 are marked in Fig. 1). F (t) = exp( At) is the 
nucleon form-factor with A = 2 GeV -2 . S 1: 5 2 are defined as 5 lj2 = 1 — k 1:2 /p lj2 ( k 1: 
k 2 , pi, p 2 , are marked in Fig. 1). The factor exp (2/3 (ti + h)) with (3=1 GeV~ 2 
[20] takes into account the effect of the momentum transfer dependence of the non- 
perturbative gluon propagator given by (p 2 is the Lorentz square of the momentum 
carried by the non-perturbative gluon): 

D(p 2 )=D e W (p 2 /^ 2 ). (3) 



The constants C and R are defined as: 

c ^^ v St£) 2 - (4) 

R = 9 J dQ 2 Q 2 exp (-3Q 2 ) = 1. (5) 

Here G and g are the non-perturbative and perturbative quark gluon couplings 
respectively, fi is the range of the non-perturbative gluon propagator (3) and D its 
magnitude at vanishing momentum transfer. From data on the elastic scattering of 
hadrons one infers D G 2 fi 30 GeV -1 and /i m 1 GeV. 

— * 

The constant R shows the structure of the loop integral. Q T is the transverse 
momentum carried by each of the three non-perturbative gluons. 

Fig. 1 seems to describe exclusive dijets production. In fact, as was clearly 
stated in the original paper on the Higgs production [2], our calculation really is a 
central inclusive one i.e. the radiation is present in the central region of the rapidity. 

2 Note the Regge- like dependence implied by Eq. (2). There are some controversies if such 
assumption is justified. The resent results [9, 10] on the inclusive DPE dijets production, however, 
are quite encouraging. 

3 It is worth to stress that formula (2) is only valid in the limit of « 1 and for small 
momentum transfer between initial and final quarks. 
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3 Differential cross-section 



Having the matrix element (2) we can write the formula for the differential cross- 
section: 

da = (2-sy 1 (2tt) ~ 8 \M pp \ 2 dPH, (6) 
where dPH is a differential phase-space factor: 

dPH = d%5 (k 2 ) d 4 k 2 5 (k 2 2 ) d A ri 5 (r 2 ) d 4 r 2 5 (r 2 2 ) 

x 6 (ki) (k° 2 ) 6 (r°i) (r 2 °) ( Pl + p 2 - h - k 2 - n - r 2 ) . (7) 

Expression (6) is to be integrated over all variables except rapidities and trans- 
verse momenta of the produced gluons. Following closely the method used in [17] 
we obtain the following result for the differential cross-section: 

In the above expression y lj2 are the rapidities of the produced gluons, u + = Ui +u 2 , 
u = (ui —u 2 )/2, C E = C 16 ^ 8 and ^1,2 = 2 ((3 + A — a'\n5i t2 ). 5i, 5 2 are expressed 
by rapidities and transverse momenta as follows: 

S iy /s= \ui\ exp (yi) + \u 2 \ exp (y 2 ) , 

S 2 Vs = \ui\ exp (-yi) + \u 2 \ exp (-y 2 ) . (9) 

The differential cross-section (8) gives a Gaussian cut-off on the total transverse 
momentum of the produced pair. Putting u + = i.e ui = —u 2 = u everywhere 
except in the exponent and performing the integration over (u + ) 2 we obtain: 

Taking into account (9) and the definition of L 12 we finally obtain: 



-2e 



^3 / A„l2 



|2 



Here Ay — yi — y 2 , y — (yi + y 2 )/2. E T — \ui\ — \u 2 \ is the transverse energy of one 
of the produced gluons. 

This completes the calculations of the differential cross-section. 
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Some comments are in order. 

(i) Since y 2 « (^#) 2 = 144 and ^cosh 2 (^) = 5i5 2 « 1 the differential 
cross section (11) depends very weekly on the sum of the dijet rapidities. It is 
worth to mention that the kinematical limit of the double pomeron exchange region 
depends on both the sum y and the difference Ay of the rapidities, as seen from (9). 

(ii) The main uncertainty in the expression (11) is the value of G 2 /4ir. Following 
[2, 16] we take it to be [21] about Lin fact it should be considered only as an order 
of magnitude estimate. 



4 Comparison with the CDF Run I results 

The CDF collaboration has presented [18] results on the central inclusive DPE dijet 
production cross sections. 

At Run I (y/s =1.8 TeV) the cross section for the central inclusive dijets of 
Ej > 7 GeV [E T > 10 GeV] is measured to be 43.6 ± 4.4(stat) ± 21.6(syst) nb 
[3.4 ± l(stat) ± 2(syst) nb]. The kinematics is following: 0.01 < 5 1 = 5 P < 0.03, 
0.035 < 5 2 = 5p < 0.095, jets are confined within —4.2 < y < 2.4 and the gap 
requirement 2.4 < y gap < 5.9 on the proton side. 

It should be noted that in the above experiment the protons were not detected 
and the DPE events were enhanced by a rapidity gap requirement on the proton 
side. Thus, in principle, there is no specified cuts on the outgoing proton 4 . 

Integrating 5 (11) over the appropriate kinematical range we obtain the results 
shown in Table 1. The running coupling constant g 2 /4:7i, appearing in (4), is evalu- 
ated at 2Ef n i.e. 0.15 and 0.14 for Ef n = 7, 10 GeV respectively. G 2 /4ir is taken 
to be 1. 



Transverse energy 


CDF 


a 


Ej > 7 GeV 


43.6 ± 26 [nb] 


70 [nb] 


Ej > 10 GeV 


3.4 ± 3 [nb] 


30 [nb] 



Table 1: Comparison of the results obtained in the present paper with the CDF 
results on the central inclusive DPE dijets production cross sections. G 2 /4ir is 
taken to be I. 

As can be seen the original Bialas-Landshoff model [2] for DPE diffractive pro- 

4 In principle the result (11) should by multiplied by a factor ^1 — exp ^ — L 1 ^^mL ) ^ ^ where 

Vgap is the maximum value of the gap. In the present case, y^f = 5.9, this factor is close to 1. 
5 Note an identical final state particle phase space factor ^ . 
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duction gives correct order of magnitude for the measured central inclusive cross 
sections. 

5 Discussion 

The results shown in Table 1 do not take gap survival effect {S 2 ap ) into account i.e. 
the probability of the gaps not to be populated by secondaries produced in the soft 
rescattering. Following [5,22] we take it for the Tevatron energy to be about 0.1 (it 
is the best we can do). Multiplying the results shown in Table 1 by S 2 gap we obtain 
the results presented in Table 2. 



Transverse energy 


CDF 


° X Sgap 

G 2 /4n = 1 


° X Sgap 

G 2 /47r = 0.4 


Ej > 7 GeV 


43.6 ± 26 [nb] 


7 [nb] 


43 [nb] 


Ej > 10 GeV 


3.4 ± 3 [nb] 


3 [nb] 


19 [nb] 



Table 2: Comparison of our results with the CDF results on the central inclusive 
DPE dijets production. Gap survival factor is taken into account. 

Some comments seems to be in order. 

(i) If we assume G 2 /An = 0.4 so that we reproduce the result for E T > 7 GeV, 
the result for E T > 10 GeV overestimates the measured cross section about 5 times. 
It means that some suppression with E T is needed. 

(ii) The expression (11) (+S 2 ap ) do not include the Sudakov factor i.e. the proba- 
bility that the rapidity gaps survive QCD radiation. We think that in our particular 
case, E T > 7, 10 GeV, the Sudakov physics is not so essential, however, the lack 
of Sudakov E T suppression can be seen. It is not obvious how the Sudakov factor 
should be included in the presented model and this problem is currently under our 
consideration. 

(iii) The factor S 2 ap is not a universal number. It depends on the initial energy 
and the particular final state. Theoretical predictions of the survival factor, S 2 ap , 
can be found in Ref . [23] . 

(iv) We work at the partonic level. Thus, in a realistic experimental situation 
our results correspond to a smaller cone where the jet axes are confined, since 
the gluon produced just on the boundary of the cone usually produces secondary 
particles outside the cone and such event is not counted. This effect would decrease 
the calculated cross section. We checked that taking the range of the rapidity as 
—3.5 < y < 1.7 (CDF —4.2 < y < 2.4) our cross-sections decrease no more than 
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10%. To cover it fully Monte-Carlo simulations are needed [14]. 

Finally, let us note that estimates in the present paper, as well as in [2, 9, 10, 15- 
17], are based on the basis of the pure forward direction. It was first mentioned in 
[24] that such approach may lead to incorrect results. 

6 Conclusions 

In conclusion, using the Landshoff-Nachtmann model of the pomeron, we have pre- 
sented the rapidity and transverse momentum distribution of central inclusive gluon 
jets produced by double pomeron exchange in pp (pp) collisions. 

We observed a distinct dependence on the difference of gluons rapidities Ay and 
the marginal dependence on their sum y. A power-law decrease approximately as 
E~ 4 - 3 , for relatively small E T production when the Sudakov physics is not essential, 
with the increasing transverse momentum is observed. We find the model to give 
reasonable predictions, at least correct order of magnitude for the measured central 
inclusive cross sections. 
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